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Phase change material (PCM) based latent heat thermal storage (LHTS) systems offer a challenging 
option to be employed as an effective energy storage and retrieval device. The performance of LHTS 
systems is limited by the poor thermal conductivity of PCMs employed. Successful large-scale utilization 
of LHTS systems thus depends on the extent to which the performance can be improved. A great deal of 
work both experimental and theoretical on different performance enhancement techniques has been 
reported in the literature. This paper reviews the implementation of those techniques in different 
configurations of LHTS systems. The influence of enhancement techniques on the thermal response of 
the PCM in terms of phase change rate and amount of latent heat stored/retrieved has been addressed as 


a main aspect. Issues related to mathematical modeling of LHTS systems employing enhancement 
techniques are also discussed. 
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1. Introduction 


The increasing gap between the global demand and supply of 
energy is becoming a major threat as well as a challenge for the 
engineering community to fulfill the needs of the energy hungry 
society. Many forums and energy management groups have been 
formed to emphasize the storage of energy in both industrial and 
domestic sectors, in any possible form. The utilization of the 
abundant source—solar thermal energy and hot waste streams 
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available in industries has attracted the scientific community to 
provide attractive solutions for the problems on energy conserva- 
tion and storage/retrieval. 

Thermal energy can be stored in the form of sensible heat in 
which the temperature of the storage material varies with the 
amount of energy stored. Water or rock can be the best example 
[1]. Alternatively, thermal energy can be stored as latent heat in 
which energy is stored when a substance changes from one phase 
to another by either melting or freezing. The temperature of the 
substance remains constant during phase change. 

Of the two, latent heat thermal energy storage technique has 
proved to be a better engineering option due to its various 
advantages like large energy storage for a given volume, uniform 
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Table 1 
Applications of LHTS units. 


S. no. Fields 


Solar based dynamic space power generation 
Solar thermal applications 

Industrial waste heat recovery 

Automobiles 

Cooling of electronics 

Textiles 

Passive heating of buildings 

Air conditioning systems 
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energy storage/supply, compactness, etc. Hence, different config- 
urations of latent heat thermal storage (LHTS) units find their wide 
applications in various engineering fields which are listed in 
Table 1. 

A summary of various practical applications of LHTS units can 
be found in Refs. [2,3]. LHTS units employ phase change materials 
(PCMs) which undergo change of phase (solid-to-liquid and vice 
versa) during the energy transfer process. During the last four 
decades many such materials, with wide range of melting/freezing 
point, have been identified and studied extensively. These 
materials include organic (e.g. paraffin), inorganic (e.g. salt 
hydrates) and salt eutectics (e.g. CaClp-MgCl2-H20). Wide range 
of PCMs with their properties, advantages and limitations have 
been comprehensively reported in Refs. [2-4]. 

Though LHTS units serve as a better energy storage device as 
they have a high thermal energy density and possess a nearly 
isothermal operation, the phase change material loaded in the unit 
possesses a very low thermal conductivity, which drastically 
affects the performance of the unit. The effect of the lower value of 
conductivity is reflected during energy retrieval or withdrawal 
with an appreciable temperature drop during the process. As a 
result, the rate of phase change process (melting/solidification of 
PCM) has not been up to the expected level and the large-scale 
utilization of LHTS units, remains unsuccessful. Therefore, to tackle 
the above-mentioned drawbacks, it becomes necessary to improve 
the thermal performance of the LHTS units employing PCMs. The 
various techniques adopted for enhancing the thermal perfor- 
mance of the LHTS units are enumerated below: 


e Using extended surfaces. 

e Employing multiple PCM’s method. 
e Thermal conductivity enhancement. 
e Micro-encapsulation of PCM. 


This paper presents a review on these performance enhance- 
ment techniques in order to provide information on relative merits 
and demerits of the various possible enhancement techniques. 


2. Using extended surfaces 


Fins or extended surfaces are used to provide additional heat 
transfer surface in thermal systems. In LHTS, the role of different 
configurations of fins in enhancing the performance has been 
studied extensively by various researchers. The LHTS, studied until 
now, are of two types: 


System packed with PCM which transfers heat with a hot/cold 
heat transfer fluid (HTF), the one used in solar water heaters, 
space based power plants, waste heat recovery systems, etc. 
System packed with PCM transfers heat with hot/cold boundary 
wall, the one used in electronic cooling. 


The location of the fins in systems where there is no HTF is 
obviously at the PCM side. But if the PCM transfers the heat with 


HTF, the location of the fins is generally based on the relative heat 
transfer coefficient. This is because, the efficiency of fin increases 
with decrease in heat transfer coefficient [5]. In most of the 
systems, the PCM side heat transfer coefficient is less than that of 
the HTF side and thus it is a general practice that fins are on PCM 
side. However, it is important to understand that whether the 
presence of fins does really give the appreciable performance 
enhancement or not. Because the governing heat transfer 
mechanisms may be different for different phase change processes 
(melting and solidification). Moreover, the heat transfer mechan- 
isms during the phase change processes depend upon the 
configuration and orientation of the system and fins. The 
proceeding discussions deal with heat transfer mechanisms and 
role of fins in determining the modes of heat transfer during 
melting and solidification separately. 


2.1. Melting 


During melting, heat is transferred to the PCM first by 
conduction and later by natural convection. This is because, the 
solid region moves away from the heat transfer surface and the 
thickness of the liquid region increases near the heat transfer 
surface. Since thermal conductivity of liquid PCM is less than that 
of solid PCM, the heat transfer by conduction almost becomes 
negligible as the melting process continues. The further melting is 
mostly by natural convection due to the density gradient that 
exists within the liquid PCM. 

The influence of natural convection on the location of melt front 
has been investigated by Lamberg and Siren [6]. The melting in a 
semi-infinite PCM storage, with an internal fin, was studied 
analytically. The analytical model consists of the well-known 
Newmann solution, which assumes that the heat is transferred 
only by conduction and natural convection can be negligible. It is 
found that though Newmann solution is exact solution, it 
underestimates the location of the melt front. 

Lamberg et al. [7] conducted numerical study on melting of 
PCM in a rectangular enclosure with and without the natural 
convection effect and the results were compared with experi- 
mental results. It was observed from the results that when the 
natural convection effect was ignored, the PCM took double the 
time of actually it takes in reality to reach the maximum 
temperature. 

Stritih [8] traced the interface location during melting of PCM in 
a rectangular storage through experimental work. The experi- 
mental results were compared with conduction based analytical 
solution. Only at the beginning of melting process, the analytical 
results showed good agreement with experimental results. Once 
melting began, the experimental interface location was much 
ahead of that one obtained from analytical calculation. 

Zhang et al. [9] experimentally studied the melting of n- 
octadecane, which was discretely heated at a constant rate from 
one side of an enclosure. During the melting process as time 
passed, the temperature in the upper region of the liquid was found 
to be higher than that in the lower region, consequently 
accelerating the melting process. It was due to the reason that 
natural convection of the liquid phase was developed, which led to 
the ascending of hot liquid and the descending of cold liquid. The 
discrete heat sources had nearly the same temperatures only at the 
early stage of the process. This indicates that the early stage of heat 
transfer is conduction controlled. As the melting process pro- 
gressed, the surface temperatures became different, increasing 
from top to bottom, which means the further heat transfer was 
convection dominated. 

Jellouli et al. [10] dealt with the melting of PCM in rectangular 
enclosure heated from below. The isotherms obtained from the 
numerical study revealed the conduction dominated melting 
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process at the early stage, as they were parallel to the hot bottom 
side of the enclosure. As time elapsed, the deformation of the 
isotherms from the horizontal became more accentuated. Since the 
deformation of the isotherms from the horizontal is a clear 
indication of effect of convection, it is concluded that the 
mechanism of heat transfer was gradually shifted to natural 
convection. 

Like in rectangular enclosure, natural convection has also been 
found as dominant heat transfer mechanism in cylindrical and 
spherical enclosures. Ng et al. [11] investigated the melting of PCM 
stored in a horizontal cylindrical annulus heated isothermally from 
inside. The natural convection was modeled by Boussinesq 
approximation. The results showed recirculation cells which are 
representations of natural convection flow in the liquid PCM. Regin 
et al. [12] analyzed the melting behavior of paraffin wax 
encapsulated in a cylindrical capsule surrounded by hot water. 
The experimental results showed the influence of natural 
convection on the melting process as the process occurred 
axisymmetrically, i.e. melting in top region was much faster than 
that in bottom region. It is also reported that better agreement 
between numerical results and experimental results could be 
achieved only after including the natural convection in the 
numerical modeling. 

Jones et al. [13] conducted experimental and numerical study 
on melting of PCM in cylinder heated from side vertical wall. 
Digital photographs, which were taken during the melting process, 
revealed that only in earlier stages of melting, the melt layer 
thickness was uniform along the vertical direction (conduction 
dominated). As the time elapsed, the melt layer thickness began to 
vary along the vertical direction. This is the indication of 
convection dominated melting. The shapes of the melt front at 
different times are similar to that of observed in rectangular 
enclosures. Akgun et al. [14] observed similar kind of heat transfer 
during the experimental work on melting of paraffin in a tube in 
shell configuration, which was oriented vertically. When the PCM 
was heated from inner wall, the molten PCM extended radially 
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outward and the molten region established a conical shape. This is 
because the convection heat transfer drives the recirculation inside 
the molten PCM. 

Though spherical system provides largest volume per unit 
surface area, only limited study could be found on their application 
in LHTS systems. However, a detailed experimental study by 
Bahrami [15] provides the details of phase change process in 
spheres. Although the solid-liquid interface was flat in the 
beginning of melting, in very short period it assumed a dome- 
like shape. The dome-like shape of interface remained unchanged 
till the end of melting. This is an indication of more melting at the 
upper portion of the spherical container. More at the upper portion 
is due to the fact that the solid sinks to the bottom as the solid PCM 
is heavier than the liquid PCM. Hence, solid PCM is in contact with 
the lower wall of the sphere through out the melting process. This 
is termed as unconstrained melting by Tan [16]. It is also reported 
that the melting at the top portion was still mainly due to natural 
convection as the natural convection cells were observed during 
the entire melting process. The experimental study by Tan [16] was 
also extended to constrained melting in which the solid PCM was 
held to prevent it from sinking down to the bottom. The 
asymmetric melting was more pronounced in constrained melting. 
The digital images revealed that the melting front took oval- 
shaped pattern during the melting due to the existence of natural 
convection cells at the top and bottom halves of the sphere. 
However, faster melting was observed at the top portion as 
stronger natural convection prevailed at top. Similarly, the melting 
process in spheres was found to be asymmetric by Ettouney et al. 
[17]. Due to temperature gradient, hotter melt at the base of the 
sphere could rise upward and loses part of heat to the cooler melt 
and solid PCM. As a result, natural convection was initiated and 
started diminishing the role of conduction. 

From these studies, it is clear that during the melting process, 
the heat transfer is dominated by natural convection. The 
domination of natural convection during melting of PCM in 
various geometries of LHTS has been illustrated in Fig. 1. As it can 


Heated wall ~~ 
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Fig. 1. Natural convection dominated melting process: (a) rectangular system [9]; (b) cylindrical system [13]; (c) cylindrical annulus [11]; (d) spherical system [16]. 
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be seen, due to natural convection more melting takes place in the 
top portion of all types of geometries. Therefore, if fins are provided 
on PCM side, then the presence of fins should not necessarily be for 
improving the conduction rate, but more importantly for improv- 
ing the natural convection process. However, when top wall of the 
rectangular system is used as heated wall, then melting would be 
dominated by conduction for longer period of time. This has been 
revealed by Lacroix and Benmadda [18] as the isotherms obtained 
were parallel to the heated top wall. The temperature gradient was 
found in the negative y-direction and thus condition became 
stable. Even the existence of buoyancy driven flows in the molten 
PCM could be noticed, they were found to be relatively weak. 
Because of the absence of natural convection, in the top wall 
heating system, the melting process took more time (approxi- 
mately double) as compared to bottom wall heating. Similarly, 
only conduction dominated melting process was observed by 
Ettouney et al. [19], when the HTF was allowed to flow downward 
in a vertical shell and tube configuration. Hence, the role of natural 
convection in the melting process can be ruled out in a rectangular 
module with top wall heating. However, if the system with top wall 
heating is coupled with environment, i.e. the bottom wall is not at 
adiabatic condition, then natural circulation in the liquid phase 
could be expected [20]. The melting of PCM in a vertical cylinder 
(since it was a two-dimensional formulation, the resulting system 
is nothing but rectangular cavity) from the top was investigated by 
Pinelli et al. [21]. Although the results of conduction based 
numerical study could show good agreement with experimental 
results, considerable difference was observed at advanced stages of 
melting. In the extended study on the same configuration by Pinelli 
and Piva [20], natural convection was accounted for. The numerical 
results showed significant quantitative improvement in the 
agreement with experimental data. This indicates that natural 
convection exists even in top wall heating system, as the system is 
coupled with environment. Reddy [22] conducted numerical study 
on the melting process in a tilted rectangular cavity, which is a 
solar integrated collector water heating system. The top wall was 
used as solar radiation receiver and through the bottom wall the 
stored heat was supplied from melted PCM to cold water. The 
simulation was performed without fins and with 4, 9 and 19 fins. 
The liquid fraction of PCM was observed after 24 h. Although all the 
finned systems showed good amount of melting as compared to 
unfinned system, only system with 9 fins was observed with 95% of 
melting. This means, the system with 9 fins provided optimum 
performance. 

Rectangular module is a basic heat storage unit that is related to 
various areas such as solar energy system, geothermal energy, 
nuclear reactors, etc., because of its simplicity [23]. Since it is 
simple in geometry and boundary conditions, extensive study has 
been made so far. In the numerical study on top wall heated 
rectangular module by Akhilesh et al. [24], the natural convection 
was ignored. As mentioned above, if the system is not coupled with 
environment, this assumption is reasonable and realistic. The 
analysis included the effect of adding more number of fins per unit 
length on the melting process. More number of fins per unit length 
led to more heat transfer area and thus energy stored. However, it 
was found that beyond a value (which is termed as critical value by 
authors) increasing number of fins could not enhance the 
performance considerably. Gharebaghi and Sezai [25] investigated 
the performance enhancement in the rectangular device due to 
horizontal fins added to the vertical heated walls. Both the heated 
walls were kept at a constant temperature and higher than the 
melting point of the PCM. The results indicated that the heat 
transfer rate could be increased due to the presence of fins. It is also 
reported that vertical heated walls with horizontal fins should be 
preferred to horizontal heated walls with vertical fins, as the 
former resulted in double the heat transfer rate as compared to the 


latter. The heat transfer enhancement factor was found to be 
increasing as the distance between the adjacent fins was decreased 
(corresponding to system with more number of fins). From the 
results, it can be interpreted that beyond some value, increasing 
the number of fins would lead to only marginal increase in heat 
transfer rate. This is because, as the number of fins is increased 
then there will be a hampering effect on the buoyancy driven flows 
and the melting becomes conduction-dominated process [26]. 

Lacroix and Benmadda [26] studied the melting process in a 
rectangular enclosure under the influence of horizontal fins 
emerging from vertical heated wall (Fig. 2). The numerical 
simulation included the effect of number of fins and length of 
fins on the melting rate. In case of shorter fins, the melting rate was 
almost independent of number fins. It was also found that once the 
melting front crossed the tip of the fins, the melting process 
became similar to that of a case without fins. Hence, the presence 
of fins was hardly felt. On the other hand, more number of longer 
fins could steadily enhance the melting rate in all cases. At higher 
values of heated wall temperature, increasing the number of fins 
produced only little enhancement in melting rate. This is due to the 
hampering effect on the buoyancy driven flows. This shows that 
the optimum number of fins depends on temperature of heated 
wall. It is concluded that few longer fins (4 fins of 0.03 m long each) 
are far more efficient in enhancing the melting rate than more 
number of shorter fins (19 fins of 0.01 m long each). Few longer fins 
could enhance the melting rate significantly even with small 
temperature difference between heated wall and melting point. 
This was found to be more efficient as compared to increasing the 
temperature of heated wall above the melting point. Similar kind of 
system was also studied by Lamberg et al. [7] and Lamberg [27]. 
Since the objective of the work was to validate the numerical 
results, no attempt was made to study the effect of number of 
horizontal fins on melting process. However, a system with two 
fins was compared with system with no fin. The experimental 
results indicated that considerable enhancement in the melting 
rate is realistic in rectangular enclosure with horizontal fins. 

In electronic devices, such as computers, digital cameras, 
cellular phones, etc., heat sinks have been in use for cooling of the 
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Fig. 2. Rectangular LHTS unit with horizontal fins emerging from vertical heated 
wall [26]. 
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Fig. 3. Fluid motion during melting [28]: (a) wider system; (b) narrow system [30]. 


components. Such heat sinks are simple rectangular modules with 
internal fins. In recent years, PCM based heat sinks (PCM is packed 
between the fins) have been studied widely. Different modules of 
PCM based heat sinks studied are heat sinks with vertical fins 
emerging from top heated surface [24], heat sinks with horizontal 
fins emerging from vertical heated surface in which portion of the 
fin tips immersed in PCM [26], heat sinks with vertical fins 
emerging from heated bottom surface [29-31], etc. These studies 
have mainly focused on performance enhancement as a result of 
adding PCM in the conventional finned heat sink rather than 
adding fins in the PCM. However, the effect of vertical fins 
emerging from heated bottom surface can still be explored from 
Shatikian et al. [30] and Shatikian et al. [31]. They studied the 
melting of PCM for different thickness of PCM stored between fins 
(wider and narrow cases) and different fin thickness. For a fixed 
size of rectangular module, wider system corresponds to system 
with less number of fins and vice versa. The results showed that the 
melting rate in the narrow cases was faster than that of in the 
wider cases. It was found that in wider cases, the fluid motion 
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became stronger and stronger as the melting progressed. On the 
other hand, in narrow cases, the fluid motion was hardly seen even 
at advanced stages of melting. The intensity of the fluid motion in 
wider and narrow cases has been illustrated in Fig. 3. 

Lacroix and Benmadda [18] also noticed that the onset of 
natural convection was gradually prevented when the distance 
between the fins was decreased. Therefore, it can be concluded that 
natural convection would exist if less number of fins are used. 
However, Lacroix and Benmadda [18] found that too large distance 
between the fins led to the reduction in total heat transfer surface 
area. This means for a fixed size of module, it is important to 
optimize the distance between the fins, i.e. number of fins. As 
mentioned earlier, Reddy [22] has found that increasing the 
number of fins from 9 to 19 caused slow melting rate. Although for 
the system studied, optimum number of fins has been declared, 
this cannot be generalized as no non-dimensional analysis was 
carried out. According to Lacroix and Benmadda [18] the optimum 
distance between the fins decreases as Rayleigh number increases 
and thus, this has been expressed as a function of Rayleigh number, 
which is valid for 2 x 10°>Ray<8.57 x 10° (Fig. 4). The 
expression suggested is 


W =aRaH+b 


where W is distance between the fins, Ray, is Rayleigh number 
(height based), a = —4.17 x 1078 and b = 1.4376. 

Shatikian et al. [30] reported that thick fins remained at the 
temperature of heated surface uniformly along the length. On the 
other hand, for the same length, thin fins experienced temperature 
gradient. From heat transfer point of view, it would be desirable if 
the fins perform at steady state. But too thick fins would lead to 
reduction in the number of fins. Therefore, the thickness of fins 
should also be optimized along with number of fins for the best 
performance of the system. 

In LHTS units employed for solar heating systems, A/C systems, 
waste heat recovery systems, space based solar power systems, 
etc., the HTF exchanges heat with PCM stored in the unit. The 
performance of different types of heat exchangers used as LHTS 
units has been investigated by many researchers. For the recent 
review on various kinds of heat exchangers proposed by 
researchers, the readers are referred to [32]. Among these shell 
and tube/concentric double pipe heat exchangers have been 
proved as high efficient for minimum volume [33]. Few studies 
have assessed the performance enhancement due to employment 
of fins in shell and tube/concentric double pipe heat exchangers. 

Lacroix [33] developed a three-dimensional numerical model 
for melting process in shell and tube LHTS with the PCM on the 
shell side and the HTF flowing inside. Annular fins were provided 
around the tube as shown in Fig. 5. Natural convection was 
included through effective thermal conductivity as a function of 
Rayleigh number in the conduction equation. The results showed 
that large amount of heat was conducted through the fins along the 
radial direction. For all ranges of mass flow rate and inlet 
temperature of HTF, considerable amount of increase in energy 
stored was observed due to the presence of fins. Increasing the 
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Fig. 4. Rectangular LHTS unit with vertical fins emerging from top heated wall [18]. 
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HTF 


Fig. 5. Shell and tube LHTS unit with annular fins [33]. 


number of fins resulted in more increase in heat transfer rate and 
thus energy stored at all conditions. But the enhancement factor 
due to fins was dependent on mass flow rate and inlet temperature 
of HTF. It is from the fact that maximum enhancement was 
observed for moderate flow rates and small inlet temperature. On 
the other hand, for a larger mass flow rate and inlet temperature of 
HTF, the enhancement factor was not very significant even with 
more number of fins. Later, similar configuration was studied by 
Zhang and Faghri [34]. The numerical study included the effect of 
fin height. It has been shown that the liquid fraction of the PCM at 
any time during the melting could be increased by increasing the 
fin height. This is because, the melting fronts on both the sides of 
the fins were greatly affected by the height of the fin although the 
effect of fins on melt front was hardly felt in between the fins. The 
investigation also focused on the effect of initial subcooling on the 
position of melt front and liquid fraction. Generally many inorganic 
PCMs show considerable amount of subcooling, which in turn 
badly affects the system performance [35]. Zhang and Faghri [34] 
too observed significant reduction in performance due to 
subcooling. The effect of subcooling on performance was almost 
zero on both sides of the fins. This clearly shows that the fins are 
very effective in countering the reduction in performance due to 
subcooling effect of PCM. 

Seeniraj et al. [36] investigated transient behavior of high 
temperature PCMs stored in finned shell and tube heat exchanger. 
Such high temperature LHTS units can be used in space based solar 
power plants. Since in solar space power systems microgravity 
conditions prevail, natural convection was assumed to be 
negligible. The study has revealed that if unfinned tube is used, 
then some quantity of PCM nearer to the exit of the tube would 
remain in solid state. This is because, nearer to the exit, the 
difference between HTF temperature and PCM’s melting point 
would be very small. A few number of annular fins are capable of 
maintaining relatively high temperature difference between HTF 
and melting point and thus melting could be found everywhere in 
the axial direction. For a fixed size of LHTS unit, significant increase 
in energy stored was observed due to the presence of fins. 
However, the increase was substantial only in the initial stages of 
melting. A decrease in percentage increase in energy stored was 
observed as the process continued. At any particular time, the 
percentage increase in energy stored was directly proportional to 
number of fins. But adding more number of fins would effectively 
affect the amount of PCM stored. Hence, it is concluded that the 
number of fins is to be suitably selected to meet the desired 
amount of energy to be stored. 

As reported by Seeniraj et al. [36], the enhancement factor due 
to fins decreases as the melting process continues. This is true only 
in case of conduction dominated melting. As mentioned earlier 
except in few cases, melting process is dominated by natural 
convection and in those cases, the enhancement factor due to fins 
would continuously increase throughout the melting process. This 
has been shown by Liu et al. [37]. They designed a new spiral 
twisted fins attached to the inner tube surrounded by a PCM stored 


outer tube. The inner tube was an electrical heating rod for 
supplying the heat to the PCM. During the pre-melting period 
which is conduction controlled, the temperature distributions in 
the PCM were more uniform than those without fins. This indicates 
faster heat transfer in case of finned system. Although the 
equivalent thermal conductivity was increased by fins during 
pre-melting period, the increase in equivalent thermal conductiv- 
ity started decreasing after some time. But once melting begun, the 
equivalent thermal conductivity started increasing once again. 
This proves that fins are very effective not only in increasing the 
melting rate, but also in reducing the pre-melting time. It is also 
shown that better enhancement is possible if fins of small width 
are used. Fine fins of same amount of material provide more 
effective area than that of large fins and thus better conduction 
rate. Moreover, fine fins divide the annular space into more 
number of smaller sections, which enhances the natural convec- 
tion. 

The foregoing discussions have focused on adding fins to the 
PCM side. If HTF of low thermal conductivity (e.g. air) is used, then 
the main thermal resistance is experienced in the HTF side rather 
than in the PCM side [38]. Therefore, Zhang and Faghri [38] 
concentrated on enhancing the forced convective heat transfer 
inside the tube by employing internal longitudinal fins. As 
compared to unfinned system, in the finned system the melting 
process was accelerated because of increased forced convective 
heat transfer coefficient and direct heat transfer from fins to PCM. 
It should be noted that the fins were more efficient only at low 
Reynolds number. The effect of dimensions and number of fins on 
the liquid fraction was also investigated. Though it was observed 
that increase in thickness, height and number of fins led to increase 
in liquid fraction at any time, beyond some point the increase was 
less significant. Therefore, it is concluded that more number of fins 
would be beneficial if thicker and longer fins are used. It is also 
reported that the adverse effect of initial subcooling of PCM can be 
countered by using internal longitudinal fins. 


2.2. Solidification 


Contrary to melting process, solidification is dominated by 
conduction. During solidification natural convection exists only in 
the beginning and as the time goes the effect of natural convection 
becomes almost zero as compared to the effect of conduction [39]. 
The heat transfer characteristics of PCM during solidification have 
been explained by many investigators. 

Stritih [8] compared the analytical results [conduction based] of 
solidification process in a rectangular storage unit with experi- 
mental results and obtained a good agreement between two. The 
natural convection was found to be 10 times lower than that in 
case of melting process. Ettouney et al. [19] studied the 
solidification in a shell and tube heat exchanger. It is reported 
that the solidification isotherms were almost parallel to the heat 
transfer surface and solidified layers took the same shape of the 
heat transfer surface. Hence, natural convection effects were found 
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to be negligible whether the direction of HTF flow was upward or 
downward. Ettouney et al. [17] obtained similar results even in 
spherical storage device. For the same module size, HTF 
temperature and HTF velocity, the heat transfer coefficient during 
melting was higher than that of during solidification. It indicates 
that melting was governed by natural convection where as it is not 
so in case of solidification. 

These studies prove that in all configurations of LHTS, during 
solidification the solidified layers are formed from heat transfer 
surface and remain parallel to it. Although natural convection 
exists in the liquid PCM at earlier stages, it diminishes rapidly as 
the solidification progresses and the process becomes conduction 
dominated. This is because as the time goes, the liquid volume 
becomes smaller and smaller. Therefore, if fins are to be employed 
to enhance the solidification rate, then the sole purpose is to 
improve upon the conduction heat transfer. 

To compare the heat transfer with and without fins in 
rectangular module, Stritih [8] estimated fin effectiveness experi- 
mentally. The effectiveness is defined as the ratio between heat 
flux with and that without fins. The results showed that the fin 
effectiveness was significantly high and because of that 40% 
reduction in solidification time was observed. Effect of vertical fins 
attached to horizontal constant temperature wall on solidification 
of high temperature PCM has been simulated by Guo and Zhang 
[40]. When the fins were not attached, solidification started on the 
wall and continuously moved away only in the vertical upward 
direction. However, after employment of vertical fins, the move- 
ment of solidified front was observed in vertical as well as in the 
horizontal direction simultaneously. More importantly, the rate of 
movement parallel to the fins was faster than that in the upward 
direction. The results also indicated that the discharge time with 
fins was almost 1/30th of that without fins. The time required for 
complete solidification was found to be decreasing almost linearly 
with number of fins. Similarly, increasing the fin thickness led to 
decrease in solidification time. There was also a critical value for fin 
thickness beyond which the effect on solidification was less 
significant. 

Shatikian et al. [41] too observed that the solidification was not 
only initiated at the horizontal wall but also at the vertical fins 
attached to the wall. It is also added that at any time heat transfer 
from the fin to the wall was much higher than that coming directly 
from the PCM to the wall. Moreover, as the time elapsed, direct 
heat transfer from PCM to wall became negligible. This shows that 
fins increase the heat transfer rate from PCM and if fins are not 
provided, solidification may be terminated after some time. As in 
the case of system with vertical fins, the solidification rate can be 
effectively increased by increasing the number of horizontal fins 
attached to the vertical wall [26]. However, the numerical study by 
Gharebaghi and Sezai [25] has revealed that heat flux is higher for 
system with horizontal fins than that of system with vertical fins 
for any number of fins. 

The study on performance enhancement due to fins during 
solidification has also been extended to cylindrical geometrics like 
shell and tube or double pipe arrangement. These studies have 
focused on both radial as well as axial fins attached to the heat 
transfer surface. Choi and Kim [42] conducted experimental study 
to determine the heat transfer characteristics of PCM during 
solidification in cylindrical storage system with radial fins. The 
results of temperature distributions showed that the temperature 
gradient at all radial positions were larger in case of finned system 
as compared to unfinned system. Thus, it is evident that there was 
an improved heat transfer between PCM and the tube due to radial 
fins. Also better axial conduction in the solid phase due to fins was 
observed, as the temperature distribution in the axial direction was 
more uniform in the finned system than those in unfinned system. 
Hence, the heat recovery was much faster due to radial fins. 


However, the effectiveness of the fins was found to be dependent 
on mass flow rate of HTF. For lower HTF mass flow rate, solidified 
layers were formed only on the small part of fins adjacent to the 
tube wall at earlier stages, though it was improved as the process 
progressed. At higher mass flow rates, solidified layers can be 
found on larger parts of fins right through the process and the best 
use of fins can be made. The effect of radial fins on enhancement 
has also been examined by Liu et al. [43]. They designed spiral 
twisted fins attached to the tube carrying HTF. It is reported that 
during solidification fins not only increase the conduction heat 
transfer, but also natural convection, which prevails at earlier 
stages. Because of this, employing fins were found to be more 
effective at earlier stages than at latter stages. Nevertheless, the 
enhancement factor was as high as 250% due to fins. It is also 
recommended to use less width fins as they produce more effective 
enhancement. 

When shell and tube or concentric double pipe heat exchanger 
is used in LHTS units, usually PCM is stored in shell/annulus and 
the inner tube surface acts as heat transfer surface. During 
solidification, the solidified PCM between the tube and liquid PCM 
increases the thermal resistance and thus poor heat transfer to the 
liquid PCM, which is away from the tube. If fins are attached to the 
tube, the liquid PCM can have a better thermal contact with the 
heat transfer surface, which in turn reduces the thermal resistance. 
The above-mentioned studies have proved this point. However, 
Velraj et al. [44] have pointed out that further enhancement is 
possible if PCM is stored in the inner tube provided with 
longitudinal fins. Fig. 6 shows the LHTS unit with internal 
longitudinal fins. This is because, in this configuration, the thermal 
resistance experienced by the liquid PCM is uniformly reduced as 
the liquid PCM farthest away from the heat transfer surface can 
have maximum thermal contact with the fins. The authors have 
also observed additional heat flow which was in circumferential 
direction. Hence due to employment of longitudinal fins, higher 
surface heat flux was observed. Moreover, the surface heat flux 
could be increased to a larger extent with more number of fins. 
Although with fins the complete solidification time was approxi- 
mately 1/n of that without fins, the effect of increasing the number 
of fins was more pronounced for larger tube radius. Recently, 
Castell et al. [45] have proved that adding vertical fins to the HTF 
side also enhances the performance of LHTS during solidification. 
This is because vertical fins on the HTF side significantly promote 
the natural convection within the HTF. The study focused on 
calculation of heat transfer coefficient through which it was found 
that the solidification time in the finned system was considerably 
less as compared to unfinned system. Therefore, it is concluded 
that in finned system same heat transfer coefficient of unfinned 
system can be achieved with less temperature difference. When 
the longer fins were used, there was decrease in heat transfer 
coefficient due to dampening effect on natural convection. 
Nevertheless, the solidification time was still less because of 
increase in heat transfer area. The LHTS unit with external 
longitudinal fins is shown in Fig. 7. 

As mentioned earlier, the number/dimensions of the fins play 
very important role in promoting the natural convection during the 
phase change process. This is very crucial in convection dominated 
phase change process. Even in conduction dominated process the 
role of number/dimensions of the fins cannot be ruled out as they 
affect the amount of material stored in the unit of fixed size. In this 
case, choice of number of fins can be a compromise between the 
enhancement factor and amount of energy to be stored. Although 
different configurations of LHTS with various fin geometries have 
been studied widely, all works have not taken into account the 
effect of number/dimensions of the fins. This can be seen from 
Table 2, which presents the list of works carried out on finned LHTS 
units. 
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Fluid bath 


Fig. 6. LHTS unit with internal longitudinal fins [44]. 


3. Employing multiple PCMs method 


Employing multiple families of PCMs in LHTS system has been 
reported as another attractive performance enhancement techni- 
que in the literature. Employing multiple PCMs means, the LHTS 
system is packed with more than one PCM of different melting 
temperatures. The heat transfer rate in LHTS unit and thus the 
performance of the system during charging (melting) and 
discharging (solidification) mainly depends on the difference 
between the HTF temperature and the melting point of PCM. If 
single PCM is used in the unit, then this temperature difference 
would obviously decrease in the flow direction of HTF. This results 
in decrease in heat transfer rate and thus poor performance of the 
unit. If multiple PCMs of different melting temperatures are packed 
in the unit in the decreasing order of their melting points, then 
nearly a constant temperature difference can be maintained during 
melting process, even though the HTF temperature decreases in 
the flow direction. This leads to almost a constant heat flux to the 
PCM. During discharging, if the HTF flow direction is reversed then 
the PCMs remain in the increasing order of their melting points. 
Once again nearly constant heat flux but from the PCM to HTF is 
possible. The use of multiple PCMs in shell and tube is illustrated in 
Fig. 8. 

The benefit of employing the multiple PCMs has been reported 
in new direction by Wang et al. [46]. The authors have introduced a 
new concept of homogeneous phase change process in which the 
phase change takes place everywhere in the PCM. This means the 
phase change rate of each component of PCM is same and remains 
constant during the process. Though a perfect homogenous PCM 


Fig. 7. LHTS unit with external longitudinal fins [45]. 


cannot be realized, it is suggested that if PCMs of different phase 
change temperatures are arranged in such a way that the phase 
change rate distribution is parabolic, then the phase change 
process would be homogeneous. It is also reported that with a 
nearly homogeneous phase change process, the phase change time 
can be reduced considerably. Hence, it can be stated that 
employing multiple PCMs seems to be a promising way to 
enhance the performance of the LHTS system during charging and 
discharging processes. 

A good amount of work has been devoted to investigate the 
performance enhancement by employing multiple PCMs in 
different configurations of LHTS units. Farid and Kanzawa [47] 
employed three PCMs of different melting points packed in 
cylindrical capsules. Air was used as HTF. During both charging and 
discharging about 10% increase in heat transfer rate was obtained 
with three PCMs as the melting/solidification started in all three 
PCMs simultaneously, where as in single PCM system the phase 
change process started at different times. Experimental study by 
Michels and Pitz-Paal [48] has revealed the benefit of using 
multiple PCMs also in shell and tube module. The shell side was 
loaded with three PCMs and synthetic oil was allowed to flow 
through the inner tube. This unit corresponds to the one used for 
solar power plant. It was observed that in case of three PCMs unit 
larger potion of PCM underwent phase change process during the 
cycle. Hence, over a period of time three PCMs unit possessed 
higher phase change fraction and storage/discharge capacity as 
compared to the unit with single PCM of higher melting point. 

In the above-mentioned configurations, HTF flows through the 
system and hence all the PCMs arranged in the flow direction (axial 
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Table 2 
Studies on finned LHTS units. 


Reference System geometry Geometry/orientation of fins Phase change Effect of no. Effect of dimensions 
process of fins of fins 
Lamberg and Siren [6] Rectangular Rectangular/between two vertical heated surfaces Melting Not included Not included 
Lamberg et al. [7] Rectangular Rectangular/between two vertical heated/cooled Melting and Not included Not included 
surfaces solidification 
Stritih [8] Rectangular Rectangular/vertical Melting Not included Not included 
Lacroix and Benmadda [18] Rectangular Rectangular/emerging from bottom heated surface Melting Included Not included 
Reddy [22] Rectangular Rectangular/emerging from top inclined heated surface Melting Included Not included 
Akhilesh et al. [24] Rectangular Rectangular/emerging from top heated surface Melting Not included Not included 
Gharebaghi and Sezai [25] Rectangular Rectangular/between two horizontal heated surfaces/ Melting and Included Not included 
between two vertical heated surfaces solidification 
Lacroix and Benmadda [26] Rectangular Rectangular/emerging from vertical heated surface Melting and Included Included (length) 
solidification 
Lamberg [27] Rectangular Rectangular/emerging from vertical heated surface Melting Included Not included 
Kandasamy et al. [29] Rectangular Rectangular/emerging from bottom heated surface Melting Not included Not included 
Shatikian et al. [30] Rectangular Rectangular, emerging from bottom heated surface Melting Included Included (thickness) 
Shatikian et al. [31] Rectangular Rectangular/emerging from bottom heated surface Melting Included Included (length) 
Lacroix [33] Shell and tube Annular/around HTF tube Melting Included Not included 
Zhang and Faghri [34] Shell and tube Annular/around HTF tube Melting Not included Included (length) 
Seeniraj et al. [36] Shell and tube Annular/around HTF tube Melting Included Not included 
Liu et al. [37] Shell and tube Spiral twisted/around HTF tube Melting Included Included 
(width, thickness) 
Zhang and Faghri [38] Shell and tube Rectangular/longitudinal, inside the HTF tube Melting Included Included 


Lamberg [39] 

Guo and Zhang [40] 
Choi and Kim [42] 
Liu et al. [43] 
Velraj et al. [44] 
Castell et al. [45] 


Rectangular 
Shell and tube 
Shell and tube 
Shell and tube 
Cylindrical tube 
Cylindrical tube 


Annular/around HTF tube 
Annular/around HTF tube 


Rectangular/between two vertical cooled surfaces 


Spiral twisted/around HTF tube 
Rectangular/internal longitudinal 
Rectangular/external longitudinal 


(length, thickness) 
Not included 
Included (thickness) 
Not included 
Included (width) 
Not included 
Included (length) 


Not included 
Included 
Not included 
Not included 
Included 
Not included 


Solidification 
Solidification 
Solidification 
Solidification 
Solidification 
Solidification 


direction) are in contact with the HTF. In solar thermal applications 
the LHTS system may be kept in the storage tank, which contains 
the HTF [49]. The PCM stored in the LHTS unit absorbs excess heat 
energy from the HTF during charging process and releases to HTF 
during discharging. In this case since the LHTS is surrounded by the 
HTF, there is no axial temperature variation within the HTF. In such 
kind of configurations, multiple PCMs should be arranged rather in 
radial direction than in axial direction in order to extract maximum 
benefit. As the centrally located PCM material in the single PCM 
unit has minimum thermal contact with the HTF, melting cannot 
be expected immediately. With multiple PCMs arrangement, it is 
possible to have lowest melting point PCM at the centre and other 
PCMs arranged in the increasing order of their melting points from 
the centre of the unit to the outer wall. This may lead to more or 
less simultaneous melting at all points as a nearly constant 
temperature difference can be maintained in the radial direction. 
One such arrangement has been investigated by Wang et al. [50]. 
They used three coaxial cylindrical tubes to store three PCMs and 
the entire unit was kept in water bath as shown in Fig. 9. The PCMs 
are arranged as explained above. The experiments were conducted 
only for charging process with different HTF temperatures. The 
results indicated that the melting time in three PCMs system was 


decreased by 37-42% as compared to that in single PCM module. 
This was due to the fact that the three PCMs unit terminated the 
melting process early by 15-25%. 

The concept of multiple PCMs has also been extended to the 
LHTS systems in which there is no HTF. As mentioned earlier, these 
systems correspond to those used for electronic cooling. In this 
case, PCMs can be arranged in decreasing order of their melting 
points from the heated wall for charging and vice versa for 
discharging. A study on one such system has been undertaken by 
Wang et al. [51]. It is reported that if single PCM is used in the 
system then the phase change temperature distribution would be 
constant. On the other hand, with multiple PCMs linear phase 
change temperature distributions are possible. In the systems with 
constant temperature boundary conditions, a proper arrangement 
of different PCMs could lead to an optimum linear phase change 
temperature distributions. This results in 25% reduction in phase 
change time as compared to that in case of single PCM unit. It was 
also found that enhancement due to multiple PCMs was dependent 
on number of PCMs. The authors have concluded that practical and 
economical number of PCMs may vary from 5 to 10. 

Multiple PCMs technique has found its application in the heat 
receivers of solar dynamic power system [52,53]. The PCM in the 


PCM 1 PCM 2 PCM 3 PCM 4 PCM 5 
Tml Tm2 Tm3 Tm4 Tm5 
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HTF flow during charging 


HTF flow during discharging 
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Tm - Melting point 


Fig. 8. Multiple PCMs in shell and tube LHTS unit. 
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Fig. 9. Multiple PCMs in LHTS unit of coaxial cylindrical tubes [50]. 


heat receiver absorbs solar thermal energy during sunlight period 
and the absorbed energy is extracted by HTF during eclipse period. 
Like in other applications, the single PCM unit has been found to be 
inefficient and does not work properly [53]. Cui et al. [53] carried 
out numerical study on the solar dynamic power heat receiver with 
three PCMs. At the end of the charging period, three PCM module 
has a maximum of 0.74 liquid fraction as against 0.68 observed in 
single PCM module. Similarly, at the end of discharging process the 
liquid fractions of three PCMs and single PCM units were found to 
be 0.08 and 0.02 respectively. This shows a significant enhance- 
ment in performance with three PCMs. Hence, it is concluded that 
to transfer the same amount of heat energy, the weight of the heat 
receiver can be reduced considerably by employing multiple PCMs 
which is very important in space. It is also stated that reducing the 
maximum temperature of the tube walls is very crucial as it 
reduces the thermal stress and improves the reliability of the tube. 
Besides, there must be minimum fluctuation in the HTF exit 
temperature in order to improve the cycle efficiency. With three 
PCMs module reduced tube maximum temperature and less 
fluctuation in the HTF temperature could be achieved. Seeniraj and 
Narasimhan [54] also observed a nearly uniform temperature at 
the HTF outlet for a longer period in the five PCM shell and tube 
arrangement. Farid and Kanzawa [47] have shown that though at 
the beginning of the phase change processes the variation of HTF 
exit temperature in the multiple PCMs unit was similar to that of 
single PCM unit, there was sharp decrease in the rate of variation in 
case of multiple PCMs unit. However, towards the end of both the 
phase change processes multiple PCMs unit provided larger 
variation in the HTF exit temperature as compared to single 
PCM unit. Gong and Mujumdar [52] have pointed out the role of 
HTF mass flow rate in the variation of HTF exit temperature. It was 


found that at lower and higher HTF mass flow rates, the 
fluctuations in the HTF exit temperature were larger in case of 
system with multiple PCMs as compared to the one with single 
PCM. This has led to the conclusion that for a fixed size of the unit, 
there must be an optimum flow rate to keep the HTF exit 
temperature fluctuation minimal. 

Although it is clear that multiple PCMs technique enhances the 
performance of LHTS systems, the enhancement factor is 
determined by the HTF inlet temperature and HTF mass flow rate 
[55]. The numerical study on the two PCMs shell and tube module 
by Adine and El-Qarnia [55] has revealed this fact. To compare the 
thermal performances of the two PCM and single PCM models 
during charging process, the term thermal storage efficiency was 
introduced. The thermal storage efficiency is defined as the ratio 
between the actual energy stored and maximum energy stored. 
With two PCMs system, a maximum storage efficiency of 0.98 was 
achieved for all HTF inlet temperatures and mass flow rates. From 
enhancement factor point of view, highest storage efficiency was 
obtained for lowest mass flow rate. When the mass flow rate was 
increased, the two PCMs system was efficient only for lowest HTF 
inlet temperature. Therefore, multiple PCMs unit is more efficient 
for low values of mass flow rate and inlet temperatures of HTF. 
Similarly the HTF inlet temperature plays important role even 
during charging. If very high difference between PCM melting 
point and HTF inlet temperature is maintained the performance 
would not be significant as higher temperature difference leads to 
excessive subcooling of PCM [52]. 

When multiple PCMs are employed to enhance the perfor- 
mance of LHTS systems, it is important to choose right combination 
of PCMs. This corresponds to appropriate difference between 
melting points of PCMs and relative proportions of the PCMs. Fang 
and Chen [56] conducted numerical study on shell and tube 
module during charging to examine the effect of different 
combinations of PCMs. Although this study is limited to three 
PCMs model, same conclusions may be extended to five or more 
number of PCMs models. Different sets of three PCMs were 
selected and in each case the difference between the melting 
temperatures of two adjacent PCMs was kept same. This difference 
between the melting temperatures was different for different 
cases. In all cases the middle PCM was same, named as benchmark 
PCM and chosen for single PCM unit. The results indicated that 
higher melt fractions and thus higher heat storage amounts were 
achieved for larger melting point differences with same bench- 
mark PCM proportion. Increase in melt fraction/energy storage was 
observed in all cases when the benchmark PCM proportion was 
decreased. This is because decrease in benchmark PCM proportion 
led to increase in proportion of PCM with lowest melting point 
which needed less time to melt. However, in all cases there was an 
optimum value for benchmark PCM proportion to achieve the 
greatest melt fraction/energy storage. The results further showed 
that the optimum benchmark PCM proportion increased with 
increase in melting point difference between the PCMs. Watanabe 
and Kanzawa [57] studied the effect of melting point difference 
between the PCMs on both charging and discharging rates. Unlike 
in the work by Fang and Chen [56], Watanabe and Kanzawa [57] 
considered the difference between the highest and lowest melting 
points of PCMs in the set of PCMs. By increasing the melting point 
difference, both melting and solidification times could be reduced. 
However, beyond a value which is stated as optimum both the 
times started increasing. It is also shown that this optimum 
melting point difference increased with increase in difference 
between HTF inlet temperature and initial temperature of the 
PCMs. They have also presented simple expressions to roughly 
calculate the optimum melting point difference in terms of number 
of transfer units (NTU). Therefore, it can be concluded that while 
employing the multiple PCMs method, selection of appropriate 
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Table 3 
Studies on multiple PCMs technique. 


Reference System geometry Phase of operation Nature of work HTF No. of PCMs combination 
PCMs 

Farid and Kanzawa [47] Cylindrical tube Charging and Numerical Air 3 All arbitrary 

loaded with discharging 

PCM capsules 
Michels and Pitz-Paal [48] Shell and tube Charging only Experimental and Oil 3 KNO3, KNO; + KCI and NaNO; (all real) 

numerical 
Wang et al. [50] Co-axial tubes Charging only Experimental Water 3} Stearic acid, Sliced paraffin and lauric acid 
(all real) 
Gong and Mujumdar [52] Double pipe Charging and Numerical He/Xe 5 LiF + CaF, (real) and all others arbitrary 
discharging 
Cui et al. [53] Double pipe Charging and Numerical He/Xe 3 LiF + CaF; (real) and all others arbitrary 
discharging 
Seeniraj and Narasimhan [54] Shell and tube Charging only Numerical Liquid Na 5 LiF + CaF,, LiF + MgF> (both real) and all 
others arbitrary 

Adine and El-Qarnia [55] Shell and tube Charging only Numerical Water 2 P-116 and n-octadecane (both real) 
Fang and Chen [56] Shell and tube Charging only Numerical Air 3 All arbitrary 
Watanabe and Kanzawa [57] Cylindrical tube Charging and Numerical Water 15/30 All arbitrary 

loaded with discharging 

PCM capsules 
Domanski and Fellah [59] Cylindrical tube Charging and Numerical Air 2 Both arbitrary 


discharging 


combination of PCMs is crucial to produce greatest performance 
enhancement. 

The above-mentioned investigations indicate the thermal 
behavior of LHTS systems with multiple PCMs under different 
operating conditions. The information provided can be very useful 
for design optimization of such LHTS systems. Recently, Verma 
et al. [58] have pointed out that thermal analysis of LHTS systems 
based on only first law of thermodynamics may not produce a 
system with highest possible thermodynamic efficiency although 
workable designs may be possible. The authors thus have stressed 
the importance of exergy analysis based on second law of 
thermodynamics. Exergy analysis has also been carried out on 
LHTS systems with multiple PCMs. 

Watanabe and Kanzawa [57] have presented guidelines to 
obtain maximum exergy efficiency from multiple PCMs systems. 
The results were obtained in terms of destroyed exergy and exergy 
efficiency at different operating conditions for both charging and 
discharging processes. For all values of HTF inlet temperature and 
initial temperature, there were optimum values of melting point 
difference between highest and lowest melting points. These 
optimum values led to lowest destroyed exergy and exergy 
efficiency. Further lowest destroyed exergy and thus highest 
exergy efficiency could be obtained only with smaller difference 
between the HTF inlet temperature and initial temperature of the 
PCMs. This is because the increase in the difference between the 
HTF inlet temperature and initial temperature of the PCMs would 
lead to increase in exergy loss to surroundings as the phase change 
process progresses. Domanski and Fellah [59] have compared the 
exergy efficiency of multiple PCMs unit with that of single PCM 
unit. The numerical results showed that the exergy efficiency of 
former was about 1.5 times that of latter. They have also proposed 
the melting point of the last PCM which is near the exit of the HTF 
should be close to ambient temperature so that minimum exergy 
loss to surroundings. 

Although employing multiple PCMs, seems to be an attractive 
technique to enhance the thermal performance of conventional 
LHTS units, using multiple PCMs of different melting temperature 
along with finned surfaces would be better enhancement 
technique as the fins provide the compartments to store different 
PCMs [36]. Most of studies on employing multiple PCMs have been 
carried out only with arbitrary chosen materials, i.e. the melting 
points of the PCMs in the combination are based on arbitrary 
choice. This shows that studies still remain on a large scale to 


investigate the performance enhancement of LHTS with real-time 
PCMs of different combinations. Table 3 summarizes the various 
works on multiple PCMs technique with the important informa- 
tion of whether arbitrary chosen or real PCMs are considered. 


4. Thermal conductivity enhancement 


Although conventional PCMs offer high energy density, slower 
rates of melting and solidification limit the potential applications 
of LHTS systems in practical systems. This is due to the fact that all 
the conventional PCMs—both organic and inorganic possess very 
low thermal conductivity ranging from 0.1 to 0.6 W/m K. Table 4 
lists thermal conductivity values along with latent heat values for 
few PCMs, which are identified as widely investigated by 
researchers so far. For comprehensive list, readers are referred 
to the review articles by Farid et al. [1], Sharma and Sagara [2] and 
Zalba et al. [3]. 

The thermal conductivity of conventional PCMs can be 
enhanced by employing high conductivity materials. The use of 
high conductivity materials to enhance the thermal conductivity of 
conventional PCMs has been dealt with in different ways by 
researchers, which can be summarized as follows: 


Table 4 
Thermal conductivity and latent heat of selected PCMs [1-3,60]. 


Name of the PCM Thermal conductivity Latent heat of 


(W/m K) fusion (kJ/kg) 
Organics 
n-Octadecane 0.35 (solid) and 0.149 (liquid) 245 
n-Docosane 022 194.6 
Paraffin wax 0.514 (solid) and 0.224 (liquid) 251 
Caprylic acid 0.149 149 
Capric acid 0.153 153 
Inorganics 
KNO; 0.5 266 
NaNO; 0.5 172 
MgCl2-6H20 0.694 (solid) and 0.57 (liquid) 168.6 
CaCl2:6H20 1.088 (solid) and 0.54 (liquid) 190.8 


Inorganic eutectics 

58.7%Mg(NO3)- 
6H20-41.3%MgCl2-6H20 

66.6%Urea-33.4%NH4Br 


0.678 (solid) and 0.51 (liquid) 132 


0.682 (solid) and 0.331 (liquid) 161 
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e Impregnation of high conductivity porous material with the 
PCM. 

e Dispersion of high conductivity particles in the PCM. 

e Placing of metal structures in the PCM. 

e Use of high conductivity, low density materials. 


4.1. Impregnation of porous material 


Porous structure may be either metal matrices made of 
aluminium, copper, etc., or naturally available porous material 
such as graphite. As far as metal matrices are concerned, recently 
Fiedler et al. [61] have found that copper matrices give 
approximately 80% increase in effective thermal conductivity as 
compared to that of aluminium matrices. More increase in 
effective thermal conductivity was found to be possible by using 
diamond coated copper matrices. 

The effect of impregnating high conductivity porous matrix of 
different thermal conductivities and different porosities with PCM 
on thermal performance of LHTS system was investigated by 
Mesalhy et al. [62]. The numerical study was carried out to 
examine the melting process in a horizontal cylindrical annulus. 
The results indicated that the performance enhancement was 
dependent on both porosity and thermal conductivity of the 
matrix. In pure PCM, the rate of melting was less at the bottom 
portion of the container as compared to that at top portion as the 
melting process was convection dominated. With high conductiv- 
ity porous matrix, the melting rate at lower portion could be 
increased. The enhancement due to porous structure was found to 
be dependent on porosity of matrix. All though for all porosity 
values higher melting rates were observed, decrease in porosity 
resulted in decrease in enhancement. Since low porosity values 
lead to higher effective thermal conductivity, theoretically there 
should be increase in performance enhancement. However, low 
porosity matrix dampens the motion of liquid PCM and thus 
natural convection within the liquid PCM. The results of this 
numerical study have proved this point as conduction like 
temperature profiles and streamlines were observed with lower 
porosity values. Krishnan et al. [63] obtained similar results from 
the numerical study on melting process in a rectangular enclosure 
filled with metal foam impregnated with PCM. In case of composite 
PCM, due to the restricted convective flow, the melt volume 
fraction at steady state was lower than that of pure PCM case. 
Nevertheless, substantial decrease in response time of the system 
was observed. Hence, it can be concluded that the porous matrix 
should possess not only high thermal conductivity but also high 
porosity for highest possible performance enhancement. 

Graphite is known for its high thermal conductivity, high 
electrical conductivity and high absorbability. The thermal 
conductivity of graphite varies from 24 to 470W/mK [64]. 
According to Haillot et al. [65], the thermal conductivity of PCM 
(0.2-1 W/m K) can be intensified by a factor of 5-100 depending 
upon the density of graphite added into the PCM. Cabeza et al. [66] 
have found that using PCM-graphite compound is the best option 
in applications like thermally stratified storage tank. The cheapest 
form of graphite is natural graphite flakes, which are composed of 
dense staked graphite layers. However, it is found that expanded 
graphite (EG) can be a better choice as it is of better absorbability 
due to high porosity. Generally EG can be prepared from natural 
graphite through chemical oxidation in the presence of concen- 
trated sulphuric/nitric acid/mixture of sulphuric and nitric acids 
and followed by drying up process in oven and finally by rapid 
heating in a furnace. When the EG is mixed with PCM, the melted 
PCM will be absorbed in the pores of the EG results in composite 
PCM. The thermal conductivity and thermal capacity of the 
composite PCM depend upon the relative proportion of graphite 
and PCM. 


Sari and Karaipekli [67] investigated thermal conductivity 
enhancement and storage enhancement due to the addition of EG 
in a paraffin based LHTS system. The PCM/EG composites were 
prepared with different mass fractions of EG. Besides analyzing the 
performance enhancement, the experimental study also focused 
on preparation of form-stable composite. The form-stable compo- 
site is nothing but composite of PCM/EG with no liquid PCM 
leakage. The composite PCM with 10% EG mass fraction was found 
to be more form-stable one as the paraffin was properly held by the 
capillary and surface tension forces of porous EG. Also the density 
of 10% EG composite was less than that of pure paraffin resulted in 
lighter LHTS unit with same volume. The thermal conductivity 
values obtained for composites with different EG mass fraction 
showed that the thermal conductivity can be increased by 
increasing the mass fraction of EG. The increase in mass fraction 
resulted in decrease in latent heat value and thus decrease in 
storage capacity. The optimum mass fraction was taken as 10% as 
this form-stable composite PCM possessed about four times higher 
thermal conductivity than that of pure paraffin with small decrease 
in latent heat value. It should be mentioned that the melting 
temperature of form stable composite was found to be lesser by 
1.4 °C. The melting timings of pure paraffin and composite PCMs 
were compared to examine the impact of improvement in the 
thermal conductivity. The results indicated significant reduction in 
melting time (about 32%) when form-stable composite PCM was 
used in the LHTS system. Contrary to the results of Sari and 
Karaipekli [67], Yin et al. [68] have reported that the thermal 
conductivity of EG/paraffin composite steadily increased only up to 
the EG mass fraction value of 6.25. Beyond this value thermal 
conductivity started decreasing. The presence of some lacunas and 
gases in the loose structure of composite PCM with high EG mass 
fraction has been cited as reason the reduction in thermal 
conductivity. This indicates that composite should be prepared 
with appropriate mass fraction of graphite in such way the lacunas 
and gases are not present in the pores. Nevertheless, with 6.25% EG 
mass fraction the composite was from-stable, which presented 
significant increase in heat transfer rate not only during charging 
but also during discharging. The effect was more pronounced 
during charging. The decreases in melting and solidification times 
were about 65.3% and 26.2% respectively as compared to those 
with pure paraffin. 

Recently, Kim and Drzal [69] used exfoliated graphite 
nanoplatelets (XGnP) to form graphite/paraffin composite. The 
XGnP, which are of less than 10 nm thick, and 15 wm average 
diameter can be obtained by pulverizing the EG using an ultrasonic 
processor. The results showed strong relation between mass 
percentage of XGnP and thermal conductivity. The continuous 
increase in thermal conductivity was observed up to the highest 
studied mass fraction of 7%. Besides, with all XGnP mass fractions, 
the latent heat values very closely matched with those of pure 
paraffin. As mentioned earlier, increase in mass fraction of 
expanded graphite would result in decrease in latent heat of 
PCM. In case of XGnP composite, since the high surface area of 
XGnP presented good dispersion in the PCM, there was no decrease 
in the latent heat. Hence, it can be stated that XGnP is more 
effective than EG as it enhances thermal conductivity with out 
reducing the storage capacity of pure PCM. 

The use of graphite to improve the thermal conductivity of PCM 
is not only limited to low melting point PCMs such as paraffin, but 
also proposed for high melting point PCMs such as salts or eutectic 
mixtures [60,70]. Generally composite PCMs are made by mixing 
the graphites with molten PCM. The same method was adopted by 
Pincemin et al. [70] to prepare eutectic salt/graphite composite. 
Because of high viscosity of salt and similarity in density between 
salt and graphite, the graphite particles could be well dispersed 
within the molten salts and resulted in an isotropic composite. 
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Fig. 10. Types of graphite. 


However, this method was found to be less economical as PCM of 
high melting point needed more energy. Also corrosion effects on 
equipment and safety issues have been reported as limitations of 
this method. To overcome this problem, the authors adopted cold 
compression method for the preparation of composite. In cold 
compression method, after mixing the solid PCM with graphite 
powder physically, the mixture was compressed at room 
temperature. Due to the compression most of the porosity was 
removed resulted in a monolithic composite. This kind of 
composite can easily be shaped as per the geometry of the storage 
system. The composite obtained from cold compression method 
possessed anisotropic characteristics because of rearrangement of 
graphite layer orthogonally to the compression axis. The thermal 
conductivity measurement showed that radial thermal conduc- 
tivity was much higher than that of axial one. For the same graphite 
quantity, the thermal conductivity of cold compressed composite 
was higher by two times as compared to that of composite 
obtained by mixing of graphite with molten PCM. Nevertheless, 
cold compressed composites were not found to be form-stable one 
as salt leakage was observed. This was due to the presence of 
impurities and mechanical stresses. 

Pincemin et al. [60] used graphite powders of different particle 
sizes to prepare the composite with eutectic salt. The thermal 
conductivity was found to be dependent not only on graphite 
amount but also on particle size. From the results obtained, it is 
suggested that for higher graphite amount, larger particle size 
should be preferred as it presents better conductive network. On 
the other hand, for lower graphite quantity, particle should be of 
smaller size. As far as storage capacity is concerned, smallest 
particle size led to larger capacity when the graphite quantity was 
more. For smaller graphite content, the effect of particle size was 
not pronounced. The results also revealed that there was a 
decrease in phase change temperature variation due to the 
addition of graphite, i.e. the composite PCM possessed almost a 
congruent melting point. This is obviously an additional benefit as 
it would lead to more or less isothermal operation. The different 
types of graphites used by the researchers to make composite 
PCMs are shown in Fig. 10. 


4.2. Dispersion of high conductivity particles in the PCM 


Using graphite made composite PCMs has been proven as 
efficient and successful way in enhancing the performance of LHTS 
systems. However, graphite composites can be prepared only 
through various mechanical/chemical processes like heat treat- 
ment, drying, mixing, grinding/compression, etc., which are time 
and energy consuming processes. Moreover, Elgafy and Lafdi [71] 
have reported that the porosity of graphite is very crucial in 
deciding the effectiveness of the composites. If graphite of small 
mean pore size is used, then there may be decrease in latent heat 
value as small mean pore size hinders the molecular motion and 
thus very difficult to impregnate the porous media with the PCM. 
On the other hand, increasing pore size reduces the capillary force 
resulting in leakage of liquid PCM. Since the high conductivity 


metal particles are free from these shortcomings, dispersion of 
high conductivity particles into the PCM is relatively simpler 
technique to enhance the thermal conductivity of PCM. 
Mettawee and Assassa [72] carried out experiments to 
investigate the enhancement in the performance of PCM based 
solar collector due to the dispersion of micro aluminium particles. 
During the charging process, the solar radiation was observed by 
the PCM stored in the container and during discharging, the cold 
water was supplied through the pipes embedded in the PCM. 
Although the heat transfer rate could be increased by increasing 
the mass fraction of aluminium, above 0.5 mass fraction, the 
increase in heat transfer rate was found to be insignificant. Due to 
increase in thermal conductivity, the charging time was decreased 
by 60% as compared to that of pure PCM. The effect was more 
pronounced during discharging, as the conduction dominated 
solidification presented more homogeneous process. To evaluate 
the overall benefit of embedding aluminium particles, the mean 
daily efficiency was calculated. The mean daily efficiency is defined 
as the ratio between heat gained by water during discharging and 
total incoming solar radiation during charging. The results showed 
significant increase in mean daily efficiency. For all mass flow rates 
of water, the highest mean daily efficiency obtained was 94% with 
composite PCM, where as with pure PCM it was only 55%. Although 
the performance enhancement due to particle addition is reported, 
the effect of adding the particles on the storage capacity has not 
been discussed in the paper. Khodadadi and Hosseinizadeh [73] 
have reported that the latent heat of PCM would decrease as the 
mass fraction of embedded particles is increased. They simulated 
the solidification of nanofluid [water + nano copper particles] in a 
square storage model. The speed of the solid/liquid interface 
during solidification was observed for pure PCM and for PCM with 
different particle mass fractions. At earlier stages there was not 
much increase in the speed of the interface due to the addition of 
particles. As the time elapsed the effect was more pronounced as 
the interface was at more advanced locations in case of higher 
mass fractions. This resulted in considerable reduction in overall 
solidification time. Besides copper and aluminium particles, silver 
has also been tried to obtain composite PCM. Zeng et al. [74] 
prepared a composite made of 1-tetradecanol and silver nano- 
particles to investigate the effect of silver particles on the thermal 
performance of PCM. It is found that the thermal conductivity of 
the PCM increased as the loading of silver particles was increased. 
Besides this, silver particles brought in additional benefit as there 
was no strong reaction between silver and organic PCM. The 
thermal stability of the composite was almost same as that of pure 
PCM. However, the phase change temperature of the composite 
was found to be slightly lesser than that of pure PCM. Although it is 
reported that latent heat decreases with increase in mass fraction 
of particles [73,74], no attempt has been made in those studies to 
investigate the outcomes of different particles mass fraction. The 
effect of particle mass fraction has been extensively studied by 
Seeniraj et al. [75]. The numerical study was carried out on shell 
and tube arrangement. The increase in particle fraction reduces the 
volume occupied by the PCM. This leads to decrease in the amount 
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Fig. 11. Types metal structures: (a) lesser rings [49]; (b) metal balls combined with metal screens [76]; (c) metal beads [77]. 


of heat stored, even though higher heat transfer rate could be 
achieved. Therefore, it is important to determine the optimum 
particle fraction, so that maximum heat transfer rate with 
minimum loss in latent heat storage capacity can be achieved. 
For the determination of optimum particle mass fraction, Seeniraj 
et al. [75] calculated ratio of cumulative energy stored with 
particles and to that without particles. The results revealed the 
existence of optimum particle fraction at which highest ratio of 
cumulative energy stored was observed. As a matter of fact, at any 
given conditions there was a range of optimum particle fraction 
and it is suggested to choose lower value in the band to ensure high 
mass of PCM and thus more amount of energy stored. The analysis 
further indicated that the optimum mass fraction for maximum 
energy stored depended upon the particle thermal conductivity 
and maximum allowable radius for interface location. For higher 
values of particle thermal conductivity, the optimum mass fraction 
was also higher. Similarly, increase in optimum mass fraction was 
observed with increase in maximum allowable radius for interface. 
The maximum allowable radius for interface is determined by the 
dimensions of LHTS unit. Hence, an appropriate combination of 
particle thermal conductivity and particle mass fraction is quiet 
important when the conductivity particles are to be dispersed in 
the PCM to enhance the performance of the unit. 


4.3. Placing of metal structures 


Placing of metal structures into the PCM has also been 
addressed as one of the thermal conductivity enhancement 
techniques by few researchers. Velraj et al. [49] placed thin walled 
hollow cylindrical steel structures into paraffin stored in a 
cylindrical LHTS system. This structure is named as lesser rings 
(Fig. 11a). The results indicated considerable reduction in 
solidification time due to the addition of lesser rings. The time 
for complete solidification with lesser rings was approximately 1/ 
9th of that without lesser rings, whereas with longitudinal fins, the 
time for complete solidification was around 1/4th of that without 
fins. One may come to a conclusion from this that lesser rings are 
more effective than fins. However, the above-mentioned reduction 
in solidification time could be achieved only with lesser rings 
occupying 20% of total volume as against fins occupying just 7% of 
total volume. This clearly shows that to get the same reduction in 
solidification time, more volume of lesser rings are required as 
compared to that of fins. It is also reported that the enhancement 
factor due to lesser rings depends upon the diameter of the 
cylindrical module, where as in case of fins, it remains almost 
constant for a fixed number of fins irrespective of diameter of the 
module. This proves that lesser rings are more effective than fins in 
larger storage systems. 

This statement is not only valid for lesser rings, but also for 
other types of metal structures according to Ettouney et al. [76]. 


In a double pipe arrangement, stainless steel balls combined 
with stainless steel screens (Fig. 11b) were placed inside the 
PCM to increase the melting rate. The results revealed the fact 
that enhancement due to metal balls cum screens purely 
dependent on diameter and number of balls. Increasing 
diameter and number of balls resulted in increase in enhance- 
ment and the enhancement factor of up to 3 could be achieved. 
As mentioned earlier, in smaller systems the enhancement 
factor may be limited. This has also been proved by Ettouney 
et al. [77]. The study was carried out to investigate the effect of 
steel balls inserted into the PCM stored in small spherical 
module (Fig. 11c). Both melting and solidification processes 
were analyzed. The results showed only 15% reduction in 
melting and solidification times due to the addition of balls. This 
is obviously much smaller than that in larger systems. 


4.4. Use of high conductivity and low density materials 


Due to relatively high density, the metal particles/metal 
structures may settle on the bottom surface of the container 
and add considerable weight to the system. Besides, according to 
Fugai et al. [78], all metal particles are not compatible with all 
PCMs. For example, with paraffin, aluminium is compatible, where 
as copper and nickel are not compatible. Similarly, aluminium and 
copper are not compatible with some salt hydrates. Hence, there 
has been a search for low-density high conductivity additives, 
which should be compatible with all PCMs. Since the densities of 
carbon fibers are relatively lower than those of metals and the 
thermal conductivities are almost equal to that of aluminium and 
copper, these can be better alternatives to enhance the thermal 
performance of LHTS systems. Moreover, carbon fibers possess 
high corrosive resistance and hence compatible with most of the 
PCMs. 

Elgafy and Lafdi [71] prepared a composite by adding carbon 
nanofibers of 100 nm average diameter and 20 um average length, 
into paraffin wax using shear mixing and melting. The study 
focused on heat transfer characteristics during solidification. The 
thermal conductivity of the composite was found to be increasing 
almost linearly with increase in mass fraction of nanofiber. As a 
result, high solidification rate was observed. Even with 1% mass 
fraction of carbon fiber, the reduction in solidification time is 
reported as 23%. There was also no reduction in the storage 
capacity due to the additives as the output power increased 
linearly with increase in fiber mass fraction. The study was 
extended to investigate the effect of surface characteristics of 
fibers on solidification rate and found higher transfer rate with 
surface treated fibers than that with fibers of untreated surface. 
The study has also revealed the importance of uniform distribution 
of fibers in the PCM to obtain further enhancement in the 
performance. 
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Fig. 12. Distribution of carbon fibers in cylindrical capsule [78]: (a) random 
distribution; (b) brush type. 


The importance of uniform distribution of carbon fibers has also 
been reported by Fugai et al. [78]. In a cylindrical capsule, carbon 
fibers were added in the PCM in two ways. In the first case fibers 
were randomly distributed, where as in the second case brush type 
fibers were used. The two cases of distribution of carbon fibers in 
cylindrical capsule can be seen in Fig. 12. The effective thermal 
conductivity with brush type was found to be three times higher 
than that with random type. This is because in brush type the fibers 
were distributed uniformly in such a way all the fibers were 
arranged in radial direction, which was the heat flow direction. For 
lower mass fraction of fibers, the randomly distributed arrange- 
ment could not present higher melting rate than that with pure 
paraffin, even though the effective thermal conductivity of former 
is greater than that of latter. This is due to the loss in natural 
convection in case of randomly distributed arrangement. On the 
other hand, the loss in natural convection could not affect the 
higher melting rate in case of brush type. This shows that the fibers 
should be arranged in such a way that they are oriented in the 
direction of heat flow. The superiority of brush type arrangement 
also in shell tube arrangement is proved by Fugai et al. [79]. They 
could achieve 30% and 20% higher solidification rate and melting 
rate respectively with 1% mass fraction carbon brush type fibers as 
compared to those with pure PCM. Hamada et al. [80] compared 
the fiber brush with fiber chips. The experiments were conducted 
for melting process and results revealed that fiber chips resulted in 
higher thermal resistance near the heat transfer tube and thus the 
over all heat transfer rate was lower than that for fiber brush type. 
According to Fugai et al. [81], when carbon fiber brushes are used 
around heat transfer tube, then there was a critical diameter for the 
brush. The authors observed that the thermal response of the 
system increased with increase in brush diameter. When the 
diameter of the brush was increased beyond a value at which the 
brushes cross the heat transfer tubes, there was no further 
improvement in thermal response. This was due to the fact that 
beyond this critical diameter, the thermal resistance between 
fibers and the tube surface was higher. Hamada et al. [82] have also 
discussed the critical diameter of carbon fiber brushes. It is stated 
that if brush diameter is too large, then the fibers leaving the heat 
transfer surface form a region where the density of fibers is low. 
This low density region would prevent the further improvement in 
heat transfer rate. Hence, it is important to make sure that there is 
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Fig. 13. Carbon fiber cloths stretched over heat transfer tube [84]. 


no low-density fiber region to get maximum benefit out of carbon 
fibers. 

From economical point of view, one may think carbon fibers are 
expensive. However, Hamada and Fugai [83] have reported that 
carbon fiber brushes could effectively reduce the cost and space. 
The experiments were carried out on carbon fiber brushes inserted 
LHTS unit used for air conditioning system. The LHTS unit without 
carbon brushes needed relatively bigger size to produce the 
required thermal output than that of LHTS with 0.8 volume fraction 
of carbon fiber brushes. On the other hand, for the same thermal 
output the cost of latter was only 75% of that of former. As 
mentioned earlier, there exists thermal contact resistance between 
carbon fiber brush and heat transfer surface. This may affect the 
heat transfer enhancement especially when the volume fraction of 
fiber is less. Very recently, Nakaso et al. [84] have proposed a new 
technique of using carbon fibers. They used carbon fiber cloths in 
place of carbon fiber brushes. The cloths made of carbon fibers 
were highly stretched over the heat transfer surface from end to 
end of the tubes. The application of carbon fiber cloth is shown in 
Fig. 13. It is reported that these cloths with high surface density can 
act as fins. The results showed that fiber cloths of 0.42% volume 
fraction presented the same thermal performance of fiber brushes 
of 0.75% volume fraction. It is clear that low volume/mass fraction 
is desirable as reduction in storage capacity due to the additives 
would be negligible. 

The review on thermal conductivity enhancement methods has 
revealed that the phase change rate (melting/solidification) can be 
increased considerably by adding high conductivity substances. In 
general, however, these additives may lead to loss of storage 
capacity of pure PCM. The loss in storage capacity limits the mass/ 
volume fraction of additives. Therefore it should be mentioned that 
while employing such additives, it is important to estimate the 
optimum mass/volume fraction of additives. Table 5 presents the 
list of works carried out to investigate the thermal response of PCM 
under the influence of various high conductivity additives. 


5. Microencapsulation of PCM 


It is also proposed that the heat transfer rate between PCM and 
source/sink can be increased by using microencapsulated PCMs. As 
the name suggests microencapsulated PCMs are microsize PCM in 
liquid/solid (core) form enveloped within a solid structure (shell/ 
wall). The shell can be made of wide range of materials including 
natural and synthetic polymers [85]. Microencapsulation of PCMs 
can be achieved by two methods namely chemical (e.g. coacerva- 
tion, complex coacervation, interfacial methods) and mechanical 
or physical methods (e.g. spray drying method). For the complete 
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Studies on PCM thermal conductivity enhancement by additives. 


Reference Type of additive PCM Phase change process Investigation on mass/volume 
fraction of additive 
Pincemin et al. [60] Natural graphite, NaNO; + KNO; Both melting and solidification Not investigated 
expanded graphite 
and ground 
expanded graphite 
Fiedler et al. [61] Aluminium/copper Paraffin Not investigated Not investigated 


Mesalhy et al. [62] 
Krishnan et al. [63] 


Haillot et al. [65] 


Cabeza et al. [66] 

Sari and Karaipekli [67] 
Yin et al. [68] 

Kim and Drzal [69] 


Mettawee and Assassa [72] 

Khodadadi and 
Hosseinizadeh [73] 

Zeng [74] 

Seeniraj et al. [75] 


Velraj et al. [49] 
Ettouney et al. [76] 
Ettouney et al. [77] 
Fugai et al. [78] 
Fugai et al. [79] 
Hamada et al. [80] 
Fugai et al. [81] 
Hamada et al. [82] 
Hamada et al. [83] 
Nakaso et al. [84] 


metal matrix 
Metal matrix 


Metal matrix 


Expanded graphite 


Graphite 

Expanded graphite 
Expanded graphite 
Exfoliated graphite 


Aluminium particles 
Copper nano-particles 


Silver nano-particles 
Metal particles 


Lesser rings 

Metal screens and beads 
Metal beads 

Carbon fibers 

Carbon fibers 

Carbon fibers 

Carbon fibers 

Carbon fibers 

Carbon fibers 

Carbon fiber cloths 


Not available 
[non-dimensional study] 
Not available 
[non-dimensional study] 
Paraffin, stearic acid, 
sodium acetate hydrate 
and barium hydroxide octa hydrate 
Sodium acetate hydrate 
n-Docasane 

Paraffin 

Paraffin 


Paraffin 
Water 


1-Tetradecanol 
Not available 
[non-dimensional study] 
Paraffin 
Paraffin 
Paraffin 
Paraffin 
Paraffin 
n-Octadecane 
n-Octadecane 
n-Octadecane 
Paraffin 
Paraffin 


Melting only 
Melting only 


Both melting 
and solidification 


Both melting and solidification 
Melting only 

Both melting and solidification 
Both melting and solidification 


Both melting and solidification 
Solidification only 


Both melting and solidification 
Melting only 


Solidification only 

Melting only 

Both melting and solidification 
Not investigated 

Both melting and solidification 
Melting only 

Melting only 

Melting only 

Both melting and solidification 
Solidification only 


Not investigated 
Not investigated 


Investigated 


Not investigated 
Investigated 
Investigated 
Investigated 


Investigated 
Investigated 


Investigated 
Not investigated 


Not investigated 
Not investigated 
Not investigated 
Not investigated 
Not investigated 
Investigated 

Investigated 

Not investigated 
Investigated 

Not investigated 


Fig. 14. Microscope profiles of microencapsulated PCMs [85]: (a) from spray drying method; (b) from coacervation method. 


procedure of these methods readers are referred to [85-87]. The 
microscope profiles of microencapsulated PCMs obtained from 
spray drying and coacervation methods can be seen from Fig. 14. 
The thermal performance of microencapsulated PCMs can be 


Table 6 


Studies on microencapsulation of PCM. 


expected to be better than that of conventional PCMs. This is 
because small PCM particles provide larger heat transfer area per 
unit volume and thus higher heat transfer rate. Besides this, 
microencapsulated PCMs bring in more advantages like less 


Reference 


Core PCM 


Shell material 


Microencapsulation method 


Comparative evaluation 


Hawlader et al. [85] 
Chen et al. [86] 
Alkan et al. [87] 
Rao et al. [94] 
Hawlader et al. [96] 
Zou et al. [97] 

Choi et al. [98] 
Ozonur et al. [99] 


Hawlader et al. [100] 


Paraffin wax 
1-Bromohexadecane 
n-Docasane 
n-Octadecane 
Paraffin wax 
Hexadecane 
Tetradecane 
Coco-fatty acid 


Paraffin wax 


Gelatin + acacia 

Amnoplastics 
Polymethylmethacrylate 
Polymethylmethacrylate 
Gelatin + acacia 
4-Diisocynate/ethylenediamine 
Melamine-formaldehyde 
Gelatin + gum arabic 


Gelatin + acacia 


Complex coacervation/spray drying 
Not available 

Emulsion polymerization 

Not available 

Complex coacervation 

Interfacial polycondensation 

In situ polymerization 
Coacervation/complex coacervation 


Coacervation 


No 

No 

No 

No 

No 

No 

No 

Yes (during both charging 
and discharging) 

No 


Table 7 


Mathematical modeling of LHTS systems employing performance enhancement techniques. 


Reference 


System geometry 


Enhancement technique 


Phase change process 


Numerical formulation 


Performance parameters investigated 


Lamberg et al. [7] 

Lacroix and Benmadda [18] 
Reddy [22] 

Gharebaghi and Sezai [25] 
Lacroix and Benmadda [26] 
Shatikian et al. [30] 
Shatikian et al. [31] 
Lacroix [33] 

Zhang and Faghri [34] 
Seeniraj et al. [36] 

Zhang and Faghri [38] 
Lamberg [39] 

Guo and Zhang [40] 


Velraj et al. [44] 


Farid and Kanzawa [47] 


Cui et al. [53] 

Seeniraj and Narasimhan [54] 
Adine and El-Qarnia [55] 
Watanabe and Kanzawa [57] 


Domanski and Fellah [59] 

Pincemin et al. [60] 

Krishnan et al. [63] 

Khodadadi and 
Hosseinizadeh [73] 

Seeniraj et al. [75] 

Fugai et al. [79] 

Hamada et al. [83] 


Nakaso et al. [84] 


Rectangular 
Rectangular 
Rectangular 
Rectangular 
Rectangular 
Rectangular 
Rectangular 
Shell and tube 
Shell and tube 
Shell and tube 
Shell and tube 
Rectangular 
Shell and tube 


Cylindrical tube 


Cylindrical tube 
loaded with 
PCM capsules 
Double pipe 
Shell and tube 
Shell and tube 
Cylindrical tube 
loaded with 
PCM capsules 
Cylindrical tube 
Rectangular 
Rectangular 
Square 


Shell and tube 
Shell and tube 
Shell and tube 


Shell and tube 


Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 
Fins 


Fins 


Multiple PCMs 


Multiple PCMs 
Multiple PCMs 
Multiple PCMs 
Multiple PCMs 


Multiple PCMs 
Additive-graphite additive 
Additive-metal matrix 
Additive-copper nano-particles 


Metal particles 
Additives-carbon fiber brushes 


Additives-carbon fiber brushes 


Additives-carbon fiber cloths 


Melting and solidification 
Melting 

Melting and solidification 
Melting and solidification 
Melting and solidification 
Melting 

Melting 

Melting 

Melting 

Melting 

Melting 

Solidification 
Solidification 


Solidification 


Melting and solidification 


Melting and solidification 
Melting 
Melting 
Melting and solidification 


Melting and solidification 
Melting and solidification 
Melting 

Solidification 


Melting 
Melting and solidification 


Melting and solidification 


Solidification 


Enthalpy/effective heat capacity 
Enthalpy 

Enthalpy 

Enthalpy 

Enthalpy 

Enthalpy 

Enthalpy 

Enthalpy 

Effective heat capacity 
Enthalpy 

Effective heat capacity 
Effective heat capacity 
Enthalpy 


Enthalpy 


Not available 


Enthalpy 
Enthalpy 
Enthalpy 
Second law based formulation 


Second law based formulation 
Effective heat capacity 
Enthalpy 

Enthalpy 


Enthalpy 
Effective heat capacity 
Effective heat capacity 


Effective heat capacity 


Phase change rate 

Melting rate and melt fraction 

Melt fraction and average HTF temperature 

Melt fraction and heat flux 

Melting rate and melt fraction 

Melt fraction and heat flux 

Melt fraction 

Melt fraction, heat transfer rate and energy stored 
Melt fraction 

HTF exit temperature and energy stored 

Melt fraction 

Melt fraction and solidification rate 

Melt fraction, heat flux and time for 

complete solidification 

Melt fraction, heat flux and time for 

complete solidification 

HTF exit temperature and heat stored/discharged 


Melt fraction and HTF exit temperature 

Energy stored and HTF exit temperature 

Melt fraction, energy stored and heat transfer rate 
Exergy efficiency and exergy destroyed 


Exergy efficiency and exergy stored 
Phase change rate 

Melt fraction 

Melt fraction 


Energy stored and effectiveness 

HTF exit temperature 

Thermal output, energy stored/retrieved and 
heat transfer rate 

Heat release rate 
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reaction of PCM with container material, ability to withstand 
volume change during phase change, etc. [85]. 

The above stated advantages in general and easy incorpora- 
tion in particular have made microencapsulated PCMs very 
attractive in improving the thermal performance of building 
materials [88-90], textile fabrics [91,92] and Thermal fluids 
[86,87,93-95]. 

Many studies have been carried out on microencapsulated 
PCMs with different combination of PCMs and shell materials like 
paraffin wax as core and gelatin+acacia as shell [85,96], 1- 
bromohexadecane as core and amnoplastics as shell [86], n- 
docasane as core and polymethylmethacrylate as shell [87], n- 
octadecane as core and polymethylmethacrylate as shell [95], 
hexadecane as core and 4-diisocynate + ethylenediamine as shell 
[97], tetradecane as core and melamine-formaldehyde as shell 
[98], coco fatty acid mixture as core and gelatin + gum Arabic as 
shell [99], etc. 

These studies mainly focused on preparation and character- 
izations of microencapsulated PCMs. Alkan et al. [87] have 
showed that microencapsules containing n-docasane have good 
thermal and chemical stability. The results of accelerated thermal 
cycle testing showed no significant change in the phase change 
temperature and enthalpy even after 5000 thermal cycles. 
Hawlader et al. [100] prepared relatively larger capsules of 
paraffin (2.833 mm diameter) and found that the storage capacity 
of capsules remained almost same even after 1000 thermal cycles. 
The geometrical structure could be retained after 1000 thermal 
cycles. Zho et al. [97] have reported that microencapsulated 
hexadecane can be suffered up to a temperature of 282°C. 
Similarly, as reported by Ozonur et al. [98], microencapsulated 
fatty acid mixture could retain the geometric profiles even after 
50 cycles and the chemical structure of the pure PCM would 
remain same after the microencapsulation. Hence, it is clear that 
the commonly used PCMs can be successfully microencapsulated 
which are proven to be thermally and chemically stable. 
However, Hawlader et al. [96] have pointed out that the ratio 
between PCM and wall material is critical in determining the 
quantity of PCM that can be encapsulated out of original PCM 
quantity. It is found that increase in wall material led to 
encapsulation of more PCM. However, this would result in less 
quantity of PCM and hence decrease in storage capacity. The 
decrease in storage capacity due to encapsulation would be more 
when physical method (spray drying method) is used as 
compared to that with chemical method (complex coacervation) 
according to Hawlader et al. [85]. 

As mentioned earlier, due to large heat transfer area per unit 
volume of the container, microencapsulated PCMs can be 
expected to present higher transfer rate than that of pure PCMs. 
Higher transfer rate results in large energy storage and release 
capacity and thus Hawlader et al. [85] and Hawlader et al. [96] 
have recommended microencapsulated paraffin for applications 
like solar energy storage system. However, the results obtained by 
Ozonur et al. [99] showed that microencapsulated fatty acid 
presented higher heat transfer rate as compared to that of pure 
fatty acid only during charging. On the other hand, during 
discharging, pure PCM was found to be with higher heat transfer 
rate. Except the study by Ozonur et al. [99], none of the other 
works on microencapsulated PCMs has compared the perfor- 
mance of the microencapsulated PCM with that of pure PCM. In 
the work by Ozonur et al. [99], the PCM used was fatty acid 
mixture which has the melting point range of 22-24 °C. Hence, the 
pure fatty acid was in liquid state at room temperature and during 
charging, the energy was stored in the form of sensible heat. On 
the other hand, microencapsulated fatty acid was found to be with 
melting point range of 29-31 °C and hence phase change could be 
observed when heated from room temperature. Therefore, this 


comparative study has not addressed the effect of microencap- 
sulation on the heat transfer enhancement during phase change 
process. Recently, Regin et al. [101] have reported that the rigidity 
of the microencapsulated PCM matrix dampers the natural 
convection and thus reduces the heat transfer rate. If a 
microencapsulated PCM should replace the relatively cheaper 
conventional PCM, then there must be a substantial heat transfer 
enhancement. As of now, it is not much clear, as none of the works 
has extensively compared the performance of microencapsulated 
PCM with that of pure PCM. The list of works on microencapsu- 
lated PCMs is given in Table 6. Hence, it can be stated that more 
work is needed especially on comparative study between pure 
PCMs and microencapsulated PCMs in order to assess the 
performance enhancement due to microencapsulation of com- 
monly used PCMs. 


6. Mathematical modeling of LHTS systems employing 
performance enhancement techniques 


Phase change problems are examples of what are closely 
referred to as moving boundary problems and their study presents 
one of the most exciting and challenging areas of current applied 
mathematical research. The existence of a moving boundary 
generally means that the problem does not admit a simple closed 
form analytical solution and accordingly much research has 
focused on approximate solution techniques. In general, phase 
change problems involve a transient, non-linear phenomenon with 
a moving liquid-solid interface whose position is unknown a priori 
and also flow problems associated with HTF. In addition, the two 
phases of PCM may have different properties and configuration of 
the LHTS unit may differ with the applications. Substantial number 
of numerical based studies on LHTS units employing enhancement 
techniques and on their performance measurement is available in 
the literature. The numerical formulations widely implemented so 
far are enthalpy method and effective heat capacity method. 

In the enthalpy method, the enthalpy which is a function of 
temperature is considered as a dependent variable along with 
temperature. Thus the enthalpy based conduction equation is valid 
for both solid and liquid phases. Besides, it is valid for solid-liquid 
interface. Therefore, there is no need to track the interface which 
makes this formulation attractive [102]. In the effective heat 
capacity method, the heat capacity of the PCM during phase change 
process (effective heat capacity) is introduced. The effective heat 
capacity is directly proportional to the stored/released energy 
during the phase change and the specific heat. Therefore, with 
effective heat capacity it is also possible to describe the non- 
isothermal phase change in the PCM. For the details of these two 
formulations, readers are referred to [103] for enthalpy method 
and [27] for effective heat capacity method. Besides these two 
formulations, Yingqiu et al. [104] have proposed another method 
called alternate thermal resistance method. 

The numerical studies on LHTS units employing performance 
enhancement technique have mainly focused on evaluation of 
melting/solidification time, heat transfer rate and amount of 
energy stored/retrieved in comparison with those of system 
without enhancement techniques. As mentioned earlier, Verma 
et al. [58] have pointed out the importance of exergy analysis of 
LHTS units. In view of that, it can be mentioned that a numerical 
study based on exergy analysis would be a better option for the 
complete evaluation of enhancement in performance under the 
influence of possible enhancement techniques. To the best of the 
authors’ knowledge very few numerical investigations based on 
exergy analysis are reported till date and this can be drawn as a gap 
in the contemporary research. Table 7 summarizes the various 
numerical investigations on LHTS systems employing performance 
enhancement techniques. 
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7. Conclusions 


The techniques employed to enhance the performance of LHTS 
units have been reviewed. Various techniques reviewed in this paper 
include using extended surfaces, employing multiple PCM’s method, 
thermal conductivity enhancement and micro-encapsulation of 
PCM. Considerable amount of work has been devoted to assess the 
effect of those techniques on thermal response of LHTS units. This 
review of those works aims at providing necessary information on 
features and limitations of each technique to enable further research 
in the area. Exergy analysis of LHTS units employing performance 
enhancement techniques may be taken as a main aspect in future 
mathematical modeling based research works. 
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